Localization of massive fermions on the baby-skyrmion branes in 6-dimensions 
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We construct brane solutions in 6-dimensional Einstein- Skyrme systems. A class of baby skyrmion 
solutions realizes warped compactification of the extra dimensions and gravity localization on the 
brane for negative bulk cosmological constant. Coupling of the fermions with the brane skyrmions 
lead to the brane localized fermions. In terms of the level crossing picture, emergence of the massive 
localized modes are observed. Nonlinear nature of the skyrmions brings richer information for the 
fermions level structure. It comprises doubly degenerate lowest plus single excited modes. The 
three generation of the fundamental fermions is associated with this distinctive structure. The mass 
hierarchy of quarks or leptons is appeared in terms of a slightly deformed baby-skyrmions with 
topological charge three. 

PACS numbers: ll.10.Kk, 11.27. +d, 11.25.Mj, 12.39.Dc 



I. INTRODUCTION 

Theories with extradimensions have been expected to 
solve the hierarchy problem and cosmological constant 
problem. Experimentally unobserved extradimensions 
indicate that the standard model particles and forces are 
confined to a 3-brane [l|, 0, H, H| • Intensive study has been 
performed for the Randall- Sundr am (RS) brane model in 
5 space-time dimensions [J Q ■ In this framework, the ex- 
ponential warp factor in the metric can generate a large 
hierarchy of scales. This model, however, requires unsta- 
ble negative tension branes and the fine-tuning between 
brane tensions and bulk cosmological constant. 

There is hope that higher dimensional brane models 
more than five could evade those problems appeared in 5- 
dimensions. In fact brane theories in 6-dimensions show a 
very distinct feature towards the fine-tuning and negative 
tension brane problems. In Refs. [ail, it was shown that 
the brane tension merely produces deficit angles in the 
bulk and hence it can take an arbitrary value without 
affecting the brane geometry. The model is based on 
the spontaneous compactification by the bulk magnetic 
flux. If the compactification manifold is a sphere, two 
branes have to be introduced with equal tensions. If it 
is a disk, no second 3-brane is needed. But still the fine- 
tuning between magnetic flux and the bulk cosmological 
constant cannot be avoided although non-static solutions 
could be free of any fine-tuning Q • 

Alternatively to the flux compactification in 6- 
dimensions, the nonlinear sigma model has been used for 
compactifications of the extra space dimensions [1, [t| ITol . 
[Tlj . As in the flux compactification, no second 3-brane 
is needed if the parameters in the sigma model and bulk 
space-time are tuned. 

Warped compactifications are also possible in 6 space- 
time dimensions in the model of topological objects such 
as defects and solitons. In this context strings [HI, Il3l . 
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EH EH, EH were investigated, showing that they can re- 
alize localization of gravity. Interestingly, if the brane is 
modeled in such a field theory language, the fine-tuning 
between bulk and brane parameters required in the case 
of delta-like branes turns to a tuning of the model pa- 
rameters fl6| . 

The Skyrme model is known to possess soliton solu- 
tions called baby skyrmions in 2-dimensional space (l7l . 
El| • In this paper we therefore consider the warped com- 
pactification of the 2-dimensional extra space by the baby 
skyrmions. We find that in the 6-dimensional Einstein- 
Skyrme systems, static solutions which realize warped 
compactification exist for negative bulk cosmological con- 
stant. Since the solution is regular except at the coni- 
cal singularity, it has only single 3-brane. Thus no fine- 
tuning between brane tensions is required. The Skyrme 
model possess a rich class of stable multi soliton solutions. 
We find various brane solutions by such multi-solitons. 

It should be noted that general considerations in the 
6-dimensional brane model with bulk scalar fields suggest 
that the mechanism of regular warped compactification 
with single positive tension brane is not possible [l9j |. 
However, the model under consideration is restricted to 
the bulk scalar field depending only on the radial coor- 
dinate in the extra space. The scalar field in the Skyrme 
model depend not only the radial coordinate but also 
the angular coordinate to exhibit nontrivial topological 
structure, which makes possible to realize regular warped 
compactification. 

Study of localization of fermions and gauge fields 
on topological defects have been extensively studied 
with co-dimension one JEE [H, HI HI and two 
P> US M M M M M\m- M &ny years ago, par- 
ticle localization on a domain wall in higher dimensional 
space time was already addressed [33j, |34| . The authors 
suggested the possibility of localized massless fermions 
on the 1-dimensional kink background in 4+1 space- 
time with Yukawa-type coupling manner. Later, local- 
ization of chiral fermions on RS scenario was discussed in 
Ref. [2(J. Analysis for the massive fermionic modes was 
done by Ringeval et.al, in Ref. [21]. For co-dimension 
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two, the localization on higher dimensional generaliza- 
tions of the RS model was studied within the coupling of 
real scalar fields [25[. Many studies have been followed 
and most of them are based on the Abelian Higgs, or 
Higgs mediated models with the chiral fermions. 

Problem of fermion mass hierarchy has been discussed 

in many articles [HIHSSHIslSBIsZi within 
the different mechanisms. They are based on the Yukawa 
coupling of the fermions and the Higgs (scalar) field. In 
Ref. [36[, the authors set up multiple branes and consid- 
ered localization of fermions on different branes in terms 
of Yukawa couplings to the Higgs field. In Ref. [13], the 
fermions have quantum numbers of the rotational mo- 
menta which are origin of the generation of fermions. 
The authors of Ref. [29( deal with this problem with some- 
what different approach. Conical singularity of the back- 
ground branes and the orbital angular momentum of the 
fermions around the branes are the key role for the gen- 
eration. In Ref. [26| hierarchy between the fermionic gen- 
erations are explained in terms of multi-winding number 
solutions of the complex scalar fields. They observed 
three chiral fermionic zero modes on a topological de- 
fect with winding number three and finite masses ap- 
pear the mixing of these zero modes. Although any 
brane localization mechanism is absent in their discus- 
sion, the idea is promising. In Ref. [13], the authors have 
taken into account more realistic standard model charges. 
Ref. [32| is about discussion of the fermion families from 
two layer warped 6-dimensions. The authors obtained 
Kaluza-Klein particles in 5 dimensions at first; those are 
finally regarded as light, standard model fermions in 4 
dimensions. 

Starting point of our approach is conventional and has 
somewhat similarity to Ref.|26j. We shall consider the 
localization of the fermions on the baby skyrmion branes 
with topological charge three. The localized modes of 
fermions are confirmed throu gh t he analysis of spectral 
flow of the one particle state [3J|. According to the In- 
dex theorem a nonzero topological charge implies the 
zero modes of the Dirac operator [3!|. The zero cross- 
ing modes are found to be the localized fermions on the 
brane. So the generation of the fermions is defined in 
terms of the topological charge of the skyrmions with a 
special quantum number called grandspin K3. There are 
different profiles of the zero crossing behavior for different 
K3, and it is the origin of the finite mass in our point of 
view. Nonlinear nature of the skyrmion fields has richer 
information than the case of Abelian string; the level 
comprises lowest doubly degenerate as well as single ex- 
cited modes, which can partially explains generation puz- 
zle of the fermions. In order to manifest more realistic 
mass structure, breaking of the rotational symmetry and 
the shape deformation of the background skyrmions is 
taken into account. 

The crucial difference of our approach from the other 
attempts is the representation of the fermions. It is based 
on our knowledge that even in the first generation of the 
fermions, they have small but finite masses. It means 



that the fermions are not pure chiral eigenstates. There- 
fore, in this article we employ the standard representation 
of the higher dimensional gamma matrices instead of the 
chiral one. 

This paper is organized as follows. In the next section 
we describe the Einstein-Skyrme system in 6-dimensions 
and derives the coupled equations for the Skyrme and 
gravitational fields. We derive a class of multi-winding 
number solutions. Some typical numerical brane solu- 
tions are shown. Formulation of the fermions in higher di- 
mensional curved space-time is discussed in Sec. III. Cou- 
pling of the fermions and the skyrmions is introduced 
in this section. Conclusion and discussion are given in 
Sec.IV. 



II. CONSTRUCTION OF THE 
BABY-SKYRMION BRANES 

A. Model 

We introduce a model of the 6-dimensional Einstein- 
Skyrme system with a bulk cosmological constant cou- 
pled to fermions [4(|. The action comprises 

S ^gravity ~t~ ^branc -t~ Sf crm ion • (I ) 

Here Suavity is the 6-dimensional Einstein-Hilbert action 

1 



S, 



gravity 



2k 2 



;R-A t 



(2) 



In the parameter n 2 — f/Mg, Mq is the 6-dimensional 
Planck mass, denoted the fundamental gravity scale, and 
A(, is the bulk cosmological constant. 

For Sbrane we use the action of baby-Skyrme model [13, 
M 



d x£br 



(3) 



with 



C 



—d M (t> ■ d M $ + — x d N $) 

+ M 2 (l+n •<£)], (4) 



brane — V Q 



where M,N run over 0, • • • ,5 and n — (0,0,1). <f> — 
(c/) 1 , (j> 2 , <f> 3 ) denotes a triplet of scalar real fields with 
the constraint 4> ■ 4> = I. The constants F,e,u are 
the Skyrme model parameters with the dimension of 
(energy) 2 , (energy) -1 , (energy) 3 , respectively. The first 
term in Eq.(j4|) is nothing but a nonlinear o— model. The 
second term is the analogue of the Skyrme fourth order 
term in the standard Skyrme model (the Skyrme model 
in 3+1 dimensions) which works as a stabilizer for ob- 
taining the soliton solution. The last term is referred to 
as a potential term which guarantee the stability of a 
baby-skyrmion. 
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f^~~ (m,») = (2,l) 




/ (m,») = <2,2) : 



FIG. 1: Typical results of the profile functions / 
(straight line), the warp metrices B,C (dashed,dotted 
line, respectively), as a function of y. 



The solutions of the model would be characterized by 
following topological charge in curved space-time 



Q= — I <f X (V,<- \ V 2 < 



(5) 



where V M means the space-time covariant derivative. Let 
us assume that the matter Skyrme fields depend only on 
the extra coordinates and impose the hedgehog ansatz 

4> = (sin f(r) cos nd, sin /(r) smn8, cos/(r)). (6) 

The function f(r) which is often called as the profile func- 
tion, has following boundary condition 



/(0) = -(to-1)7t, lim/(r)=7r 



(7) 



where (to, n) is arbitrary integer. This ansatz ensures the 
topological charge 



Q = n(l - (-l) ro )/2. 



(8) 



We consider the maximally symmetric metric with van- 
ishing 4D cosmological constant, 



ds 2 = B 2 {r)r ]i , v dx l "dx y + dr 2 + C 2 {r)d9 2 



(9) 



where t]^ is the Minkowski metric with the signature 
(— , +, +, +) in our convention and < r < oo and < 
9 < 2tt. This ansatz has been proved to realize warped 
compactification of the extra dimension in models where 
branes are represented by global defects [4l| . 

SWmion is the action for fermions coupled with the 
skyrmions and the gravity; we shall describe it in Sec lIIII 

The general forms of the coupled system of Einstein 
equations and the equation of motion of the Skyrme 



model are 



Gain — k 2 (— Ab<?MN + Tmn) ■ 

-^d N (y/=dF 2 $x d N $ 
J— a \ 



(10) 



+^g—d M <j)(d M <j> • (0 x d N <j>)) +/i^xn = 



(11) 



where the stress-energy tensor Tmn is given by 



Tmn = 



MN 



9mnC> 



branc 



= F 2 d M $- d N $+ ^g AB (d A $xd M $) ■ (d B $x d N $) 

+9MN £branc ■ (12) 

Inserting Eq. ([6]) into Eq. (J4J) , one obtains the Lagrangian 

n 2 sin 2 / 



•Cbrane = — B^CF^e 2 



c 2 



+ 2/1(1 + cos/) 



(13) 



where we have introduced the dimensionless quantities 
Xf, = eFx^, y = eFr, C = eFC, p, = —=5 A* (14) 



F 2 



and 



u = 1 + 



n sin / 
C 2 



(15) 



The prime denotes derivative with respect to the radial 
component y of the two extra space. The Skyrme field 
equation is thus 



1 

2u 



n 2 sin 2/ 
C 2 



where 



u C 2 + n 2 sin 2 / 



(l + /' 2 ) + 2/2 2 sin/ 



& 

/'sin 2/- 2-^ sin 2 / 

c 



(16) 



(17) 



Within this ansatz, the components of the stress-energy 
tensor (fl"2")) becomes 



T^ = -F*e 2 B 2 Vfil ,T (y) , 



m 2 7^ sin / 
r (y)=-f + 2d2 



+ /2 2 (l + cos/) (18) 



T rr = -F e T r (y) , 



it 
2' 



9 • 2 r 

n sm / _ 2 



Tr {y) = -^/' 2 + - ;r J +/2 2 (l + cos/)(19) 



Tee = ~F*C z Te(y) 



a 



2C 2 



_ / , • // 2 M • cos/) (20) 



FIG. 2: Typical results of the stress-energy tensors To,r r ,To 
(dashed, dotted, and dot-dashed line, respectively) and the 
topological charge density q (straight line). 



FIG. 3: Typical results of the curvature invariants 
R, RabR ab , RabcdR abcd and CabcdC abcd 
(straight,dashed, dotted and dot-dashed line, respectively). 



where 



u=l- 



n 2 sin / 
CP 



(21) 



The Einstein equations with bulk cosmological constant 
are written down in the following form 

3b' + 6b 2 + 3bc + c +c 2 = -a(A b + r (y)) (22) 
66 2 + 46c = -a(A b + T r (y)) (23) 
4&' + 106 2 = -a(A b + Tg(y)) (24) 

where a = k 2 F 2 is a dimensionless coupling constant 
and A;, = A b /e 2 F A is a dimensionless bulk cosmological 
constant. Also, we introduce b = B'/B, c = C'/C for 
convenience. 



B. Boundary conditions 

At infinity, all components of the energy-momentum 
tensor vanishes and the Einstein equations (f2"2"|) - (f2"4")l are 
then reduced to 

36' + 6b 2 + 3b c + c! + c 2 = -aA b (25) 
6b 2 + 4bc= -aA b (26) 
46' + 106 2 = -a~A b . (27) 

The general solution has been obtained in Refs. [HI, [l6| 
which is given by 



AeiPy-e-%P y „ bp 2 3 f 

b = V = = , C= ; —0 

Ae^y + e-iPy 2b 2 



(28) 



where A is an arbitrary constant and 



P 



-aA b 

"To - 



(29) 



Since we are interested in regular solutions with warped 
compactification of the extra-space, the functions B and 
C must converge at infinity. This is achieved only when 
A b < and A = with the solution 



B 



ei e 



C 



£2e 



(30) 



where e\ and e 2 are arbitrary constants. Then, the 
asymptotic form of the metric which realizes warped com- 
pactification is given by 



ds 2 = e x e 



•q^dx^dx 11 



dy 2 



e 2 e 



J d9 2 



(31) 



The 4-dimensional reduced Planck mass M p i is derived 
by the coefficient of the 4 dimensional Ricci scalar, which 
can be calculated inserting the metric ([9]) into the ac- 
tion ©, 



^ J d 4 xy/^)R^ = ^ J d e x^B- 2 (r)R^ 
= if f d4a: ^/-9^ i^(43 / drd6B 2 (r)C(r) 



drB 2 (r)C(r) / d^xyf-g^R^ 



2ttMI 



where the superscript (4) represents a tensor defined on 
the 4 dimensional submanifold. Thus, we find the rela- 
tion between M v i and Mq as 



M I = 2nM* 



drB 2 (r)C{r). 



(32) 
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The requirement of gravity localization is equivalent to 
the finiteness of the 4-dimensional Planck mass. For the 
asymptotic solution (|3ip. the localization is attained. 

Let us consider the asymptotic solutions for skyrmions. 
we can write 



f(y) = f + Sf(y) 



(33) 



where for y ^> 1, / ~ 0. The linearized field equations 
are given by 



6f"-5 P 6f-fiSf = 0. 



(34) 



Assuming that / falls off exponentially, one obtains for 



Sf(y) -I, with q= V^+^-s P (35) 

where f c is an arbitrary constant. 

Following regularity of the geometry at the center of 
the defect is imposed 



B'(0) = 0, C(0) = 0, C"(0) = 1 



(36) 



and we can arbitrarily fix B(0) = 1. Boundary conditions 
for the warp factors and the profile function at the origin 
are determined by expanding them around the origin. 
For the different topological sectors, the first few terms 
are schematically written down as 

f(y) = -(m - 1)tt + / (n) (0)y" + 0(y n+1 ) (37) 
b(y) = By + 0(y 3 ) (38) 
C(y)=y + Cy 3 + 0(y 5 ) (39) 



where 



(m,n) = (1, 1) 

6=-|(A h + 2A-i/'(0) 4 ), 

C = ^(A b + 2M-2/'(0) 2 -^/'(0) 4 ) (40) 
(m,n) - (1,2) 

B=--(A 6 + 2m), C = -(A fc + 2A) (41) 
(m,n) = (2,1) 



s = -f(A b -if'(o) 4 



C = --(A b -2/'(0) 2 -^/'(0) 4 



a 
12 

(m,n) = (2,2) 
B 



ol - a ~ 

4 Ab ' C= 12 Ab 



(42) 



(43) 



Thus one finds that the only /'(0) or /"(0) is the free 
parameter vicinity of the origin. 

Consider linear combinations of Eqs. (|22| -(|24 | . we ob- 
tain 

1 - a 

6' + 46 2 + 6c = --oA b + -("n r + r ) , (44) 
46c + c' + c 2 = -ioA b + ^(4r + r r - 3r e ) . (45) 



Integrating Eqs. (|44"|) . (l45p from zero to y c , we get 
S 3 (»c)B'(l/c)C'(i/ c ) 



a r 
= ~2 Ab 



B 4 Cdy - -(fir + fig) . 



(46) 



B\y c )C\y c ) 



£? 4 (7dy- | (4^o + Mr -Me)- (47) 



(l4"6"| is the 6-dimensional analogue of the relation deter- 
mining the Tolman mass whereas Eq. (|47p is the general- 
ization of the relation giving the angular deficit. Com- 
bining these the following relations are obtained in the 

y c —> oo 



Jo C 



+ f' 2 )dy = 1 



(48) 



or 



a I B 
'o 



oo rin 2 
4 



n 2 sin 2 / 



C 



2A b C + 2p,C(l - cos /) dy=l. 

(49) 



These conditions are used for checking the numerical ac- 
curacy of our calculations. 

In order to study the singularity structure of the bulk 
solutions, the several curvature invariants are computed 
[l5j |. The explicit form for the metric ((9]) are given in 
Ref.[!3] and they are 

R := 20b 2 + 86' + 2c' + 2c 2 + 86c 
RabR '■= 

806 4 + 206' 2 + 2c' 2 + 2c 4 + 646'6 2 + 4c 2 c' + 286 2 c 2 
+ 32 6 3 c + 866'c 



b'c' + 8b' c 2 



2 c + 8bcc + 86c 3 



RabcdR 



ABCD 



4c 4 + 406 4 + 166 2 c 2 + 86^6' + 4c 2 + 326^6' + 166' 



2 -/ 



I P.'2 



C A 



BCD 



C 



ABCD 



^(b'-c' + bc-c 2 ) 2 . 
o 



C. Numerical analysis 



The equations ([I5 ]) .(f2"2 ]) -([2"g ]) should be solved numer- 
ically since they are highly nonlinear. The simple tech- 
nique to solve the Einstein-Skyrme equations is the 
shooting method combined with the 4th order Runge- 
Kutta forward integration [43[. However, a unique set 
of boundary conditions at y = produces 2 distinct so- 
lutions, one of which grows exponentially and another 
decays exponentially as y — > oo. This causes instability 
of solutions when the forward integration is performed. 
Instead, we employ a backward integration which is used 
in Ref. 15j where the 6-dimensional vortex-like regular 
brane solutions were constructed. The backward integra- 
tion method requires a set of boundary conditions at in- 
finity. We, however, truncate and take the distance y ma x 
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FIG. 4: Parameter space for typical solutions which exhibit 
gravity localization around the skyrmions with Q = 3. 
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FIG. 6: Fermion number density for the background 
skyrmion with (m, n) = (1,1) with the model parameters: 
a = 0.5, A b = -0.1, p. = 0.220525915 and /'(0) = 
0.7292182131078184. The results for the coupling constant 
M — (decoupled) and M — 0.86 are plotted. If the fermions 
couple to the skyrmions, only the state K3 — is localized 
on the brane core. The other states are not observed because 
they are strongly delocalized. 

III. FERMIONS 



FIG. 5: The effective potential Vu defined by Eq.(|A2)l with 
the parameter, w = 1.1,/ = for (m,n) — (1, 1) is shown. The 
figure clearly exhibits Volcano shape. 



far enough from the origin so that the Skyrme profile 
would fall off before it reaches y max - The set of boundary 
conditions at y — y max produces a unique solution which 
satisfies the boundary conditions at y = and hence it is 
numerically stable. We present our typical numerical re- 
sults in Figs. |l|2l3] , Fig. [4] shows the fine-tuning surface 
in the model parameter space (a, ft, A;,) corresponding to 
the gravity localization condition. We should stress that 
once the desired solutions are obtained, the parameters 
are no longer arbitrary, i.e., a constraint h(a, fl, Ab) = 
is emerged. 



A. Basic formalism 

The action of the fermions coupled with the Skyrme 
field in a Yukawa coupling manner can be written as 



crmion 



with 



fen 



d 3?£fermion 



^{iT A D A - Mt- 



(50) 



(51) 



The 6-dimensional gamma matrices T A are defined with 
the help the vielbein e A and those of the flat-space 7°, i.e., 
T A := e^7 a . The covariant derivative is defined as 



D 







(52) 



where 



ie aB V Ae b B are the spin connection with 



generators := j[~fa,%]- The simbol d implies that 
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FIG. 7: Spectral flow of the fermion energy of (m, n) = (1, 1). 
The zero crossing spectrum (bold line) has the quantum num- 
ber Ks = 0, which corresponds to the localizing zero mode. 
Exchange between a heavy (dotted line) and a light level 
(straight line) , which has been discussed in Ref.[46l] in some- 
what different context, is also found. 



tp d <f) = tpd(p — (dip)4>. Here A, B = 0, • • • , 5 are the 
6-dimensional space-time index and a, b = 0, • • • , 5 cor- 
responds to the flat tangent 6-dimcnsional Minkowski 
space. The vielbein is defined through qab = e^egs — 
Vab e A e B- We introduce the form of vielbein which was 



used, e.g., at Ref.[31| 



el=B{r)8l, y, = 0, 



.3, 



4 = C(r) sin^, e% = C(r) cos 6 . 



(53) 



Straightforward calculation shows that the nonvanishing 
components of the corresponding spin connections are 



< = 5*B'{r) cos 9, Lof = 5»B'{r) sin 0, 
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= 1-C"(r) 



(54) 



The Dirac equation is 



1 9 F> f 1 C 1 ' 

+ *(cos0 7 4 + sin0 7 5 )(d r + — - -^-) 



1 



— i(sin#7 — cos6 l 7 )—dg — MT 



* = 0. 



(55) 



right or the left. The massive modes emerge as their mix- 
ing states. Our approach, however, is somewhat different 
for treatment of the masses of the fermions. Actually, 
even in the first generation the fermions have intrinsic, 
finite masses. So, treating the massive fermionic modes 
directly, we employ the standard representation of the 
higher dimensional gamma matrices instead of the chiral 
one. The eigenvalues of the Dirac hamiltonian are esti- 
mated for background brane solutions with large parame- 
ter space, and the zero modes appear as the zero crossing 
points. The standard representation of the gamma ma- 
trices in 6-dimensional can be defined as 



r 



r 



7 





7A 


—Ha 



T ■ 

—Hi 




h 









,7 



h 

-h 
h o 



a 
-a 



(56) 



where /„ means the identity matrix of the dimension 
n. The 6-dimensional spinor W can be decomposed into 
the 4 dimensional and the extra space-time components 
^(x M ,r, 9) = ^{x^) ® U(r,8). Here the 4-dimensional 
part ^(x 11 ) is the solution of the corresponding Dirac 
equation on the brane 



(57) 



in which the eigenvalues w correspond to the 4- 
dimensional masses of the fermions. The Dirac equa- 
tion in 6-dimensions is transformed as the 2-dimensional 
eigenproblem with the eigenvalue w. The following re- 
placement of the eigenfuncations greatly simplifies the 
equation of motion [28| 



u(y,9) := 

2\n B(y) 



exp 



1 



dy' 
C{y') 



The eigenproblem becomes 

H2U = wu 
where the hamiltonian is 



U(y,6). 



(58) 



Ho :=B 



Mf-4> 

'a + f ) 



e (d y 

-Mr 



ide \ 



(59) 



Here we have introduced the dimensionless coupling con- 
stant and the eigenvalue M := M/eF and uj := oj/eF. 
One easily confirms that H2 commutes with "grandspin" 



A, := h 



2_ 

2 



(60) 



The Dirac gamma matrices should satisfy the anti- where j _ f a { th ut } m momentum in the 

extra space and for convenience, we introduce 7 s := I ® 
a 3 . Thus the eigenstates are specified by the magnitude 



commutation relations {7^, 7 s } = 2?] AB and there are 
the possible candidates preserving such Clifford algebra. 
In most of previous studies in the 6-dimcnsions, they are 
based on the localization of the chiral fermions on the 
Abelian vortex. In the chiral representation, the spinors 
are expanded into the right- and the left-handed compo- 
nents and the zero mode appears as a eigenstate of the 



of the grandspin, 



A 3 = 0,±l,±2,±3-- - , 

a-3=444-, 



for odd n 
for even n . 



(61) 
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FIG. 8: The typical brane solution with (771,71) = (1,3). The parameters are a = 0.11, A b = -0.1, ju = 0.49588, / (3) (0) = 
0.1913003760462098. 



If we consider "time dependent" Dirac equation 

id T u(T, r, 6) = H 2 u(t, r, 0), (62) 

the equation is invariant under "time-reversal" transfor- 
mation T, i.e., t — > — r, 



TH 2 T- 1 =H 2 , T:=i-f 5 ®T 2 C, 



(63) 



Here C is the charge conjugation operator. Since the 
hamiltonian is invariant under time reverse, the states of 
±K% are degenerate in energy. 

Once the desired eigenfunctions are obtained, angular 
averaged fermion densities on the brane can be estimated 
as follows 



(64) 
(65) 



(p) := / CMrfU = N(y)p(y) 
p(y) := / d0u\y,0)u(y,6) 



where 



N(y) = cxp 



-41nB(y) + J" 



dy' 



C(y>) 



(66) 



By using the asymptotics of the warp factors (|30)l . we eas- 
ily find out the behaviour N(y) at zero and infinity [28j ] 



N(y ~^0)^y 
N(y^w)^- 



exp(4py + —e pv ) . (67) 

P£2 



At far from the core, the densities of the excited modes 
are more enhanced than the localized mode. 



The eigenequation (|58|) can be recast into a 
(Schrodinger like) second order differential equation. We 
assume a form of eigenfuntion 



u(y,0) := 



(68) 



/ m(yy w \ 

r]2{y)e l{l+n)e 

\ Uyy (l+n+1)9 ) 

where I is an arbitrary integer. In order to eliminate £1 
and £2 from Eq. (|58|) . we write 



(69) 
(70) 



£1 = [3(G-D l _ Vl - M sin fD l + n m ) 
6 = -P{M sin fDim - G+D l + n m ) 



13:= 
G± : 



M 



72 



where 

, -2 , -1 

(71) 
(72) 
(73) 



\B 2 



B 



D± ■■=d y ±-~ 



± M cos / 



and then, we obtain 
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/3G-dj + (5 L G- + (iG'_)d v +IG--S 



w+i 



G- 



G'_C-G-C 



C 



-(^-"M sin / + 0M cos ff')d y - ( M sin / - 5< +1 M sin / 



C 2 
l + n 
C 



I | S>rn +<j-/3Msin/c\ 2 



cos /f'C- sin /C" -.. 
(3 11 — = 1 — M(l 

c 



m = o 



i i ic' 

(S n+1 M sin f + f3M cos //'H + ( M sin/ - <5 i+ ™ +1 M sin - f3M cos ff'^r + @M sin / — 



(iG+d 2 + (5 1 G+ + PG' + )d y + (G+ - 5 l+n+1 G+ 



I + n 
C 



-(3G 

r 



. l + n ^ (l + n)C\ . 



(74) 
(75) 



where J' := (3' + -~(3. Furthermore the following replace- 
ment simplifies Eqs. ([74]) ,(|75 |l 



C(y) ■= exp 
x(y) ■■= exp 



n-Ax 



K-A 2 ■ 
7, dy 



dy m(y), 
m{y), 



(76) 
(77) 



where TZ := 
the form 



C ' 



We finally obtain the equation of 



-d 2 

V 



Kd y + Vu -A 2 d v + V 12 

AlOy + V 2 1 -dl - Tidy + V 22 



X 



= . (78) 



The explicit forms of Ai, Ai, Vij = 1,2) are given 
in appendix [A] Here we show the typical example of 
the effective potential Vn in FigfS] The shape of the 
potential is like famous Volcano type. So we expect the 
existence of localized mode of the fermions. 



B. The numerical method 

Instead of solving the second order differential equation 
(|75)) . we study the eigenproblem ([55]) directly. We employ 
the method which was originally proposed by Kahana- 
Ripka [38| for solving the Dirac equation with non-linear 
chiral background. According to the Rayleigh-Ritz vari- 
ational method (44j , the upper bound of the spectrum 
can be obtained from the secular equation; 



det (A - eB) = 



(79) 



with 



Aij — 



J cPxiplHipj, Bij = J (Pxtplvj 



where {ipi} (i = 1, • • • , N) is some complete set of the 
plane- wave spinor basis. For N — > 00, the spectrum e 
becomes exact. Eq. (|T9")) can be solved numerically. For 
simplicity, we construct a plane- wave basis in large circu- 
lar box with radius D as a set of eigenstates of the flat, 
unperturbed (B = 1, B' = 0, C = r, f = tt) Hamiltonian 



7 b 7 4 d 4 " 7 7 <% - TMt a 



(80) 




FIG. 9: Fermion number density for the background (m, n) = 
(1,3) skyrmion with the model parameters: a — 0.11, At, = 
—0.35, jl — 0.957648. The case of the coupling constant M = 
1.0. The states K3 — 0, ±1 are localized on the brane core. 
Doubly degenerate K3 = ±1 states are more localized (they 
exhibit the non-zero value at the origin). 



The solutions are 

C P =M ki 



Jp-l 



(hy)e 



U Z. J 9+l 



l(q+l)6 



(81) 



with 



M ks = 



Here oj^ 
M) and e k 



e ki + M 
e h + M 



J p -i(kiD) 

J q +i(hD 
,+ 



1-1/2 



-1/2 



-sgn(e k )k/(e k 
±\/ k 2 + M 2 . The momenta (i 
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FIG. 10: Spectral flow of the fermion energy for the back- 
ground (m,n) = (1,3) skyrmion with the model parameters: 
a = 0.11, A b = -0.35, p, = 0.957648. Two zero crossing 
spectra have the quantum number K$ = (straight line), 
and Kz = ±1 (dotted line), respectively. 
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FIG. 11: Fermion number densities for the background 
(m,n) = (1,3) skyrmion about four regions depicted in 
FigHD] 



are discretized by the boundary conditions 



J p (hD) = Q, J q (kD) = 



(82) 



i 71 1 71 

where p := K 3 + - - -, q := K 3 - - + -. The orthogo- 
nality of the basis is then satisfied by 



D rD 

drrJ v {kir)J v {kjr) = / drrJ L ,±i(k t r)J„± 1 (k : jr) 



D 



1 



= &ij-;r[J v ±i{k i D)Y, v:=Ks±-T» 



(83) 



Expanding the eigenstates of Eq. ([55|) in terms of the 
plane-wave basis, the eigenproblem reduces to the sym- 
metric matrix diagonalization problem. A special care 
is taken for the estimation of the matrix element of the 
hamiltoman (|59| . In order to hold the Hermiticity of the 
matrix, the following differential rule is imposed 



\d u 



l 



dydBCiy)- ^d y ct>-(d y ^)ct> . (84) 



Size of the radius D is chosen so as to wrap the whole 
branes. Apparently the typical value D = 10.0 is suf- 
ficient because one easily observe that the brane profile 
functions, the stress energy tensors, and the curvature 
invariants are flat at y > 8 (see Figs. |ll2l3] ). 

Fig[5] shows the densities p(y) for the background 
skyrmion of (m,n) — (1,1). We display a tower of the 
massive modes together with the ground state. As one 
easily observe that only the lowest mode is peaked on the 
brane while all other modes escape from the core; there- 
fore the massive modes cannot be observed on the brane 
core. In Fig[7l we plot the fermion spectra as varying 
the coupling constant M. Only the lowest isolated mode 
decays from positive continuum to the negative. We will 



confirm that the spectrum corresponds to the brane lo- 
calized mode. 

Libanov,Troitsky have discussed relation between the 
topological charge and the fermionic generation [2(|. 4- 
dimensional fermions appear as zero-modes trapped in 
the core of the global vortex with winding number three. 
We shall investigate this speculation with our solution 
for (m,n) = (1,3). The results of the expansion at the 
origin ([57 ]) - ([55] ) are 

B = -~(A 6 + 2/2), C = ^(A b + 2p) (85) 

and f( 3 '(Q) is the shooting parameter. The typical ex- 
ample of solution of the brane is shown in FigJSl As 
is expected, we observe three localized solutions (FiglH]) 
which can be regarded as the generations of the fermions. 
A main difference from the previous studies is that we use 
non-linear type of soliton solutions for constructing the 
branes. They have richer information of the topology and 
the spectra exhibit doubly degenerate ground states and 
a higher excited state. On the other hand, the linear type 
of solitons like Abelian vortex have the spectra with dou- 
ble degeneracy only. The mass spectra of the fermions 
in our universe comprises doubly degenerate plus single 
level with very high energy; it suggests the underlying 
topology that we employs. 



C. Fermion level crossing picture 

In FigJTl we show that only one isolated mode dives 
from positive energy to negative. This behavior is called 
as spectral flow or level crossing picture (38|. Spec- 
tral flow is defined as the number of eigenvalues of 
Dirac Hamiltoman that cross zero from below minus the 
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FIG. 12: The brane solutions with large parameter space are plotted. 



number of eigenvalues that cross zero from above for 
varying the properties of the background fields. Ac- 
cording to the Index theorem a nonzero topological 
charge implies the zero modes of the Dirac operator 
[39l ] . The number of flow coincides with the topologi- 
cal charge and the zero modes are emerged when they 
cross the zero. The level crossing picture extensively 
studied in the Dirac equation with non-linear chiral back- 
ground [H, |45| , with the Higgs field of the Abelian- 
Higgs model |4rl H3] and with the non-trivial gauge fields 
(e.g.,instanton,meron) [48|, H^, [5(| • In a electroweak the- 
ory, one level crossing of the fermion in the background 
of the spphaleron barrier is observed (5lj . Sometimes 
the mechanism can be thought of as a quantum mechan- 
ical description of fermion creation/annihilation. Inter- 
estingly, we can observe the mixing of some levels, i.e., 
the energy levels cannot cross and the excited particle 
changes with the light one (see Fig(7]) . This behavior in- 
dicates that these fermions interact with each other via 
somepotential. This has been thoroughly discussed in 
Ref . [4q with somewhat different models. In this subsec- 
tion, we shall examine the level crossing behavior of the 
fermion in the brane-skyrmion background. 

Studying the spectral flow argument, the authors al- 
ways investigate evolution of the spectrum by one pa- 
rameter, e.g., time (which characterize the size or the 
strength of the background field). Our model, however, 
has many parameters which define the basic property 
of the branes. Therefore, we explore the fermion level 
behavior for the brane solutions with large parameter 
space. FigfTOl shows the result of the spectral flow of the 
solutions with (m,n) — (1,3) as a function of the cou- 
pling constant M. Fig llll is the corresponding fermion 
densities. One easily observes that the density exhibits 
various localization behaviors in every domain. In or- 
der to examine more thorough parameter dependence of 
the spectral flow, first we prepare several varieties of the 
brane solutions with large parameter space (see Fig fT^)) . 
Once the brane solutions are obtained, the model pa- 
rameters (a,p.,Ab) are no longer arbitrary, i.e., a con- 
straint h(a,jl,A},) = exists. So if we fix A(,,M, either 
a or /2 is identified as a evolution parameter. We dis- 
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FIG. 13: Spectral flow of the fermion energy of (m, n) = (1, 3) 
for M = 0.72, Afc = —0.2 and the different values of p, (a is 
the function of /I, i.e., ot(p,)). Every dot corresponds to the 
different background brane solutions. Especially, the bullet 
(•) is the lowest, doubly degenerate states and the square 
(■) is the single mode, respectively. 



play in Figs. |13lH] . the spectral flow corresponding to 
these parameters. Fig |15l shows more general spectral 
flow "cascade", i.e., the flow as functions of (p,,M). 



D. Mass splitting of the generations 

Because of the time-reversal symmetry, the spectra 
contain doubly degenerate states and, the first two gen- 
erations of the fundamental fermions should be observed 
as the degenerate states within our framework. Some 
effects can split the degenerate states. For example, in 
order to manifest the symmetry breaking, we introduce 




FIG. 14: Spectral flow of the fermion energy of (m, n) — 
(1,3) for M = 0.72, A{, = —0.2 and the different values of a. 
The bullet (•) is the lowest, doubly degenerate states and the 
square (■) is the single mode. 

the explicit mass term for the fermions of the form 

H m := B 7 § m, m := f™ 1 ° J (86) 

H m commutes with the grandspin operator and then the 
eigenstate of the hamiltonian with this term still are la- 
beled by K3. On the other hand, it breaks the time- 
reversal symmetry, thus, this additional term successfully 
splits the degeneracy. Without loss of generality, we can 
set mi = 0, m 2 = Am, which is treated as a free param- 
eter. 

Another trial is based on a fine structure of the back- 
ground solitons. In Ref.[12], the authors extensively 
studied the structure of the multisolitons with charge 
1 < Q < 5 in the baby Skyrme model. They use one pa- 
rameter family of the potentials U — /i 2 (l — <fo) s - For the 
Q = 3, the soliton has spherical symmetry below some 
critical value s (our calculation is the case of s = 1). The 
symmetry is broken above the critical value, and the so- 
lution exhibits only Z2 symmetry. So, we expect that our 
branes slightly deform and the degenerate spectra split. 

For the 3+1 Skyrme model, construction of the 
ansatz for the multisolitons with spherical coordinates 
(r, 9, <j>) called Rational map ansatz has been proposed 
by Houghton, Manton and Sutcliffe 53]. Although the 
solitons exhibit the complex platonic symmetries, they 
become tractable by using this ansatz. Coupling of the 
fermions to the Rational map skyrmions was already 
studied in Refs.[45l.l54|. In the baby Skyrme model, how- 
ever, no analytical form of ansatz for such non-continuous 
symmetry exists, then we treat the effects as a small per- 
turbation. 

In general, for the solution with non spherical symme- 
try, the profile function is modified f(r) — > F(r,ff). In 



FIG. 15: Spectral flow "cascade" of one of the zero crossing 
mode in (m, n) = (1,3), for fixed Af, = —0.2. 

the case of Z 2 , if the deformation is small, F(r, 9) can be 
expanded 

F(r, e)~/(r) + |A {r)(e 2W +e-™) (87) 

where the first term of the right hand side f(r) corre- 
sponds to the profile function with spherical symmetry. 
One can easily confirm that the solution has Z 2 symme- 
try, i.e., F(r,9 + tt) = F(r,9). Substituting into Eq.flBJ), 
the term f ■ <f> — > f ■ <p + V e thus hamiltonian (|59p can be 
written as 

H 2 -» H 2 (e) = H 2 + B^MV t (88) 
where the potential V e is defined as 
V £ (r,9;e) := 

/ - sin/(r)(e 2i9 + e- 2lB ) cos/(r)( e - ie + e^ 6 ) \ 
€ \ cos/(r)(e ie +e 5i9 ) sin f{r)(e 2W + e - 2W ) ) ' 

(89) 

for the charge three [n — 3). Since we have no de- 
tailed information about the deformation, we suppose 
e := as a constant parameter describing measure 

of the deformation. The potential breaks both the time- 
reversal and the grand spin symmetry, then the coupling 
between states with different grandspin Kz occur. Mix- 
ing of states (K 3 , K' 3 ) with large AK 3 := \K 3 — K' 3 \ can 
be negligible for small deformation. Here we compute the 
eigenproblem H2{e)u = urn, taking into account only the 
coupling with AK3 ^ 2. 

At present status of our model, we admit that it is a toy 
model for understanding the level structure of the realis- 
tic standard model fermions because we have no explicit 
realistic charges. We are able to fit just one quark/lepton 
sector in terms of a set of model parameters. First we 
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evaluate the mass of the "down sector" of quarks (d, s, b). 
FigJTS] shows the spectra. We employ a brane solution 
with the parameters a = 0.08, A 6 = -0.21, Ji = 1.4029, 
and the coupling constant M — 0.94. For obtaining the 
dimensionful values, we tentatively choose the parame- 
ters of the Skyrme model, Fe = 27800 MeV. We obtain 
masses of the down, the strange, and the bottom quark 
m d = 5.1 MeV, m s = 98 MeV, m b = 4200 MeV, for 
e = 0.012. For the lepton sector, we have to employ a dif- 
ferent solution which is characterized by the set of param- 
eters a = 0.07, A b = -0.21, jj, = 1.8838 and M = 1.03. 
The Skyrme model parameter Fe — 10506 MeV and 
when the deformation parameter is e — 0.03383, we 
obtain the masses m e = 0.5 MeV, to m = 107 MeV, 
rn T = 1779 MeV (Fig lT7|) . These are in quite good agree- 
ment with the corresponding experimental observations. 

Here, it is worth mentioning the physical implications 
of our parameter choice. In Secjlll we have computed 
the brane solutions with wide range of the parameters 
(a, Af,,fM), which are physically equivalent. The coupling 
constant M has been chosen in order that the lowest en- 
ergy levels cross zero. The splitting of the first two gener- 
ations and the third can be adjusted in terms of the choice 
of (a, Ab,p,). Essentially, the deformation parameter e is 
not a free parameter; it should be determined uniquely in 
terms of a consistent calculation of the Einstein-Skyrme 
systems. Thus, one can say that the splitting of the first 
two generations is also controlled by the brane parameter 
choice. Finally, the skyrmion parameter Fe ~ 10 4 MeV 
guarantees the size of the branes as r ~ 10" 1 fm, which 
is consistent with our observation; it is sufficiently small 
so as not to observe any evidence of the extra dimension. 



IV. CONCLUSION 

In this article, we have proposed new brane solutions 
in 6-dimensional space-time and have discussed about lo- 
calization of the fermions on them. The brane have con- 
structed by a baby-skyrmion, a generalization of a 0(3) 
non-linear a model, and have realized the regular warped 
compactification in 6-dimensional anti-de Sitter space- 
time. The metric is non-factorizable with the warp fac- 
tors either exponentially diverging or decaying in static 
solutions. But only exponentially decaying warp factors 
can allow gravity localization near the brane in the sense 
that the 6-dimensional Planck mass takes a finite value. 
Such solutions could be obtained numerically if we im- 
pose suitable boundary conditions at the distance far 
from the brane and integrate in backward. On the other 
hand, the forward integration method is unstable for any 
boundary conditions at the origin. 

After the solutions of the skyrmion branes were suc- 
cessfully found, we have studied the fermion localization 
on them. The Dirac equation in 6-dimensional curved 
space-time has been constructed in terms of introduc- 
ing the vielbein and 6-dimensional generalization of the 
gamma matrices. To treat the massive modes, we have 
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FIG. 16: Mass splitting of the down series of quarks (d, s, b) 
in terms of effect of the deformation of the background baby 
skyrmions as function of the parameter e which determines 
the strength of the deformation (in MeV). The experimental 
values are quoted from Ref . [55l] . 



used the standard representation of the gamma matri- 
ces. In order to solve the eigenproblem, we have intro- 
duced the plane wave basis in a large circular box with 
radius D. Studying identification of the charge of the 
background skyrmions and the generation number of the 
standard model fermions in our universe, we have inves- 
tigated the skrymions of the topological charge Q = 3. 
This conjecture has been confirmed through the discus- 
sion of the fermion spectral flow. It embodies the local- 
ized modes where the number equals to the background 
topological charge. As a result, we have found localized 
fermion modes on the skyrmion branes corresponding to 
their topological charge. For (to, n) = (1, 1), we have ob- 
served the solutions with some parameter ranges, which 
means the existence of the massive modes as well as the 
massless one. For (m,ri) = (1,3), three solutions local- 
ized on the brane have been found. They comprise dou- 
bly degenerate lowest modes of plus single excited state. 
This level structure well describes the experimental mea- 
surements of fermion masses. On the other hand, in the 
Higgs mediated models, the spectra exhibit only double 
degeneracy. 

As was found in Ref. (52| . the minimal energy config- 
uration of the baby skyrmion with charge three has no 
spherical symmetry, rather, it exhibits Z2 symmetry. The 
small shape deformation has an effect on the lowest de- 
generate state and splits them. We have successfully ob- 
tained the tower of the "down" sector of quark generation 
within one parameter family of the deformation. We have 
treated the deformation of the skyrmion in a perturba- 
tive way and have neglected the effect for the gravity. In 
order to address the deformation of the brane solutions 
from spherical geometry, we need to examine full Simula- 
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FIG. 17: Mass splitting of the leptons (e, \x, r) in terms of 
effect of the deformation of the background baby skyrmions 
as function of the parameter e which determines the strength 
of the deformation (in MeV). The experimental values are 
also from Ref . 15511 . 



tion of the Einstein-Skyrme system without imposing any 
symmetrical ansatz for the geometry. We believe that the 
solution will clarify the origin of the level structure of the 
standard model fermions. 

In Ref.[2l[, the authors had an elaborate analysis 
for the massive modes in a 5-dimensional anti-de Sitter 
space-time. They set a domain that is defined by the 
effective potential, and then solved the Dirac equation 
with inside and out. By imposing the boundary condi- 
tions at the surface, they obtained several massive modes. 
In our case, we also introduce a domain particularly for 
the numerical reason. We have confirmed that even if 
the radius of the domain goes to its infinity, there is no 



notable effect to the zero mode. 

Our model is a simple toy model for understanding the 
level structure of the realistic standard model fermions. 
We observed generation of masses of the down sector 
(d, s, b) or the lepton sector (e,/i,r) with good agree- 
ment of experimental observations. Perhaps most serious 
drawback of our approach is that we have much number 
of free parameters and need the independent parameter 
sets to fit each of the quark/lepton sectors. However, 
this may be justified if one considers that the quarks and 
leptons are localized on different branes but share the 
4-dimcnsional spacetime. Another defect of our model 
is of course lack of the explicit spin and charge of the 
fermions. Taking into account the effects of the break- 
ing of the flavor symmetry and assigning the realistic 
charges for all quark/lepton sectors to our solutions will 
be a great advance toward full understanding of the gen- 
eration mechanism of our universe. 

Our choice of the representation in this article might 
need more thorough consideration. Although there are 
several advantages to treat the massive modes, mecha- 
nism of localization of left-chiral fermions on the brane is 
not clear. Estimation of the asymmetry of the chirality 
on the brane is undoubtedly important subject, at least, 
about the lowest zero crossing level. The result imple- 
menting this will be reported in the forthcoming paper. 
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APPENDIX A: THE EFFECTIVE POTENTIAL 



The explicit form of Vij In Eq. (|78|) are given by 



V n = 

Vl2 = 

V n = 

V 22 = 



-/T 1 + (V +1 /r 1 - P^M 2 sin 2 / + /3G + G'_ - \(3M* sin 2//') I - !°L 
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-PM 1 -^ sin / (g + cot //' - G' + |g + ) - A 2 ^-^ 
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